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Reproducibility of a noisy limit-cycle oscillator driven by a random piecewise constant 
signal is analyzed. By reducing the model to random phase maps, it is shown that the 
reproducibility of the limit cycle generally improves when the phase maps are monotonically 
increasing. 

§1. Introduction 

When a spiking neuron receives a randomly fluctuating input, its reproducibility 
of spike generation improves compared with the case of a constant input. This phe- 
nomenon can be interpreted as phase synchronization between uncoupled nonlinear 
oscillators that receive a common fluctuating input, because repeated measurements 
on a single oscillator using the same input is equivalent to a single measurement 
on an ensemble of uncoupled identical oscillators. In our previous studies, we ana- 
lyzed the cases where the fluctuating input is given by a random telegraphic signaP^ 
or by a random impulsive signal. '^^ In this proceeding, we analyze the case where 
the fluctuating input is a slowly varying, piecewise constant random signal using 
the phase reduction technique,'*^' as a generalization towards a full treatment of 
realistic continuous random signals. 

§2. Fluctuation-induced phase synchronization 

We consider an ensemble of N identical uncoupled limit-cycle oscillators subject 
to a common fluctuating input: 

X,{t) = G{X,{t))+I{t) (2-1) 

for i = 1, • • • , A^, where Xi{t) represents the internal state of the i-th oscillator at 
time t, G{X) the intrinsic dynamics of each oscillator, and I{t) a fluctuating input 
common to all the oscillators. The fluctuating input I{t) is a piecewise constant ran- 
dom signal that takes one of M values Im G {-^i, • • • ,Im} with equal probability. 
The changes of I{t) occur at time {ti,t2, ■ ■ ■} following a Poisson process of mean 
interval r. We assume r to be sufficiently larger than the period of the oscillator. 
At each t„, I{t) changes its value in a stepwise manner. Namely, if I{tn — 0) = Im, 
its new value I{tn + 0) after the change is either of Im+i or Im-i with equal prob- 
ability. The probability density function (PDF) of the interval T„ = tn+i — tn 
between changes obeys an exponential distribution P{T) = exp(— T/r)/r. For 
each value of Im-, Ea. ()2-l|) is assumed to have a stable limit-cycle solution, whose 
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basin of attraction is the entire phase space except some unstable fixed points. 

Though our theory itself is a general 
one, we use the FitzHugh-Nagumo model as 
an example, where X = {u,v}, G{X) = 
{e(f + a — f — — n}, and I{t) = 
{0,I{t)}. The parameters are fixed at a = 0.7, 
b = 0.8, and e = 0.08. I{t) takes one of M = 7 
values € {0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2}. We 
set the mean interval between changes at r = 40 
and consider N = 25 oscillators. In the numer- 
ical simulation, small Gaussian-white noise of 
zero-mean and intensity D = 10~^ is indepen- 
dently applied to each variables of the oscilla- 
tors to incorporate the effect of external distur- 
bances. Figure H^a) displays a typical realiza- 
tion of the piecewise-constant signal. Fig. E^b) 

zero-crossing events of the v-component from v < U to v > U under tJie constant m- 
put, and Fig. ^c) zero-crossing events under the fluctuating input, sufficiently after 
initial transients. Due to the independent Gaussian-white noises, the zero-crossing 
events occur randomly under the constant input as shown in Fig. E^b) , whereas phase 
synchronization induced by fluctuating input can clearly be seen in Fig. ^c). 

§3. Reduction to random phase maps 

The phase synchronization is the result of the stabilization of each limit-cycle 
oscillator against phase disturbances due to the fluctuating input. To analyze its 
mechanism, we reduce our model to random phase maps. We consider the single- 
oscillator problem, because the stability is a property of individual oscillators. 

Corresponding to the M values of I{t), the orbit of our model moves among M 
limit cycles. Since r is assumed to be large, the orbit is on one of those limit cycles 
most of the time, except for short transients between limit cycles after the changes 
of the input, as shown in Fig. |2 Following the standard procedure, ^^'^^ we define 
a phase variable 6m{X) £ [0, 1] using the limit cycle m corresponding to the input 
Im for each m = 1, • • • , M, where and 1 represent the same phase. We specify the 
value of I{t) by m hereafter. When the input is m, i.e., I{t) = Im, the dynamics of 
the orbit can simply be described as 9m{t) = uJm by using the corresponding phase 
variable 6m, where ojm is the angular velocity of the limit cycle m. 

When the input changes from m to m', the orbit of our model originally at 
phase 9m on the limit cycle m will be mapped to new phase 9m' on the limit cycle 
m' . We describe this mapping by 9m' = Fm-*m'{Gm)-, which we call a "phase map". 
It is a periodic function on [0,1] satisfying Fm-^m'{9m + 1) = Fm^m'{9m) + 1 = 
Fm^m'{Gm), where and 1 should be interpreted as the same phase. Figure 01 
displays the phase maps of the FitzHugh-Nagumo model obtained for all contigu- 
ous pairs of (m,m'). The curves are appropriately shifted to adjust their origins. 



(a) 




Fig. 1. (a) Typical realization of a 
piecewise-constant random signal, 
(b) Zero-crossing events under a 
constant input, (c) Zero-crossing 
events under a fluctuating input. 
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We use the time step n rather than the real 
time t in the following discussion, which is the 
number of changes in I{t) from the beginning. 
Since we consider a Poisson process, the time 
step n roughly corresponds to the real time t as 
n ~ t/r, because the mean inter-impulse inter- 
val is r. Let us represent the temporal sequence 
of I{t) by m{n), and consider a situation where 
the input changes from m(n) to a new value 
m(n + 1) at t = t„ and keeps this value until 
t = tn+i for an interval of T„ = tn+i — t^- The 
corresponding dynamics of the orbit from phase 
^m(n)(™) on the limit cycle m(n) to the new 
phase ^m(n+i)(^ + 1) on the limit cycle m{n + l) 
can be described using the phase map F as 




Fig. 2. Typical trajectory of the model 
(bold arrow). M limit-cycle orbits 
corresponding to the M input values 
are shown in the background. 



^m(n+l)("' + 1) — ^m(n+l)Tn + Fm(n)^m{n+l){(^m{n){f^)) i 



(3-1) 



where ^m{n+i)Tn represents constant increase of the phase on the limit cycle m{n+l). 
Since T„ is a random variable, this equation describes random phase maps. 

§4. Stability against phase disturbances 

The stability against phase disturbances can be characterized by the average 
Lyapunov exponent of the random phase maps, Eq. H3-l() . Let us consider a small 
phase deviation Z\0^(-„)(n) from 6'rn.(n)('^)- Its linearized evolution equation is 

^0^(^n+l){n + 1) = -?4(„)^„^(„+l)(^m{n) (4-1) 

i{&m) = dFm^m' {(^m) / dOm- Therefore, the phase deviation grows as 



where 

|^^m(n)(n)/Z^e,r.(0)(0) 



n-1 



n K^m{n')-m(n'+l)(^m(n')(^')) =i exp (An) , (4-2) 



n'=0 



where we defined the average Lyapunov exponent as A = (log \ F^_^^,{6m)\) ■ The av- 
erage should be taken over all possibilities of (m, m') and over the phase distributions 
on all limit cycles. 

When the mean interval r is sufficiently large, the phase distribution on each 
limit cycle tends to be uniform, because the jumps between the limit cycles occur 
irrespectively of where the orbit is, leading to complete randomization of the phase. 
Under this condition, we can make a general statement on the sufficient condition 
for the phase synchronization: when all phase maps Fm^m'i^m) (^re monotonically 
increasing non-identity functions, the Lyapunov exponent A is negative, leading to 
fluctuation-induced phase synchronization. 
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Actually, when F!^^^,{6m) > holds for 
all m, we can bound the Lyapunov exponent A 
from above as 
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Fig. 3. Phase maps Fm-tm' {S) between 
two successive values of the input 
values m and m' . The inset is an 
enlargement showing the correspon- 
dence between the curve and the 
pair of input values (m,m'). 



where the summation is taken over all combina- 
tions of m and m', and #(m, m') represents the 
number of them. In the above inequalities, we 
utilized the fact that log F' < F' — \, and that 

Pm~*m'iX) -^m— >m' (0) — 1 becaUSe F^^fyi' (^diT^) 

is a phase map. The equality holds only when F^^^,{9m) = 1 for all m, namely, 
when the phase maps are trivial identity maps. For the FitzHugh-Nagumo model 
with the parameter values assumed here, all the phase maps -?*m^m'(^m) '^^^ mono- 
tonically increasing as can immediately be seen from Fig. Therefore, by applying 
a piecewise-constant random signal with large mean interval r, fluctuation-induced 
synchronization occurs as demonstrated in Fig. I2tc). In general, as long as the 
separation between neighboring values of I{t) are small, the phase maps should be 
monotonic, and fluctuation-induced synchronization should occur. 



§5. Summary 



We analyzed fluctuation-induced phase synchronization among uncoupled noisy 
oscillators for the case of a slowly varying, piecewise-constant random input. By 
reducing the model to random phase maps, we gave a general sufficient condition for 
the phase synchronization. Extension of our current analysis to a realistic continuous 
random signal will be tackled in the future. 
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